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Abstract 

The two-loop electroweak corrections to the anomalous magnetic moment of the 
muon, generated by fermionic loops, are calculated. An interesting role of the top 
quark in the anomaly cancellation is observed. New corrections, including terms of 
order G^amf, are computed and a class of diagrams previously thought to vanish 
are found to be important. The total fermionic correction is —(23 ±3) x 10 -11 which 
decreases the electroweak effects on g — 2, predicted from one- loop calculations, by 
12%. We give an updated theoretical prediction for g — 2 of the muon. 
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The anomalous magnetic moment of the muon, a M = (g^ — 2)/2, provides a 
precision test of the standard model and potential window to "new physics" effects. 
The current experimental average [|lj 

a e ^ = 116592300(840) x 10~ n (1) 

is in good accord with theory and constrains physics beyond the standard model 
such as 0, [| supersymmetry, excited leptons, compositeness, etc. A new experiment 
|4j] being prepared at Brookhaven National Laboratory is expected to reduce the 
uncertainty in a^ xp to below ±40 x 10~ n , more than a factor of 20 improvement. At 
that level, electroweak loop corrections become important and new physics at the 
multi-TeV scale is explored. 

To fully exploit the anticipated experimental improvement, the standard model 
theoretical prediction for should be known with comparable precision. That re- 
quires the confluence of calculational effort involving very high order QED loops, 
hadronic loop contributions, and even two loop electroweak effects. Indeed, the con- 
tributions to a M are traditionally divided up into 

a„ = a Q ED + a» adronic + af (2) 

The QED loops have been computed to very high order H 

oed = « + 0.76585738(6) (-Y + 24.0454(4) f-' 

ZTY \7T/ \7T 



+126.14(43) (-) + 930(170) f-J (3) 

where in the calculation of the r lepton loops we used m T = 1777 MeV. Employing 
a- 1 = 137.0359895(61) recommended by the PDG gives 

a QED = n 6 5 8 4708(5) x lO -11 (4) 

The uncertainty could be further reduced by a factor of 2, if we chose to use a as 
determined from the electron g e — 2, a -1 = 137.03599222(94); however, in either case 
it is well within the ±40 x 10 -11 experimental goal. 

A recent reexamination of hadronic vacuum polarization at the 0(a/n) 2 level, 
utilizing e + e~ — > hadrons data via a dispersion relation, gives 

aJ adronic (vac. pol.) = 7023.5(152.6) x 10~ n (5) 

Unfortunately, the uncertainty has not yet reached the hoped for level of precision. 
However, it is anticipated that ongoing improvements in e + e~ — ► hadrons data 
near the p meson resonance (which weighs heavily in (|Sp) and theoretical input in 
the higher energy region will significantly reduce the uncertainty. Nevertheless, the 
goal of going below ±40 x 10 -11 remains a formidable challenge. 

The result in (f|) must be supplemented by higher order, 0(a/n) 3 , hadronic vac- 
uum polarization effects 0, |3| 

a Hadronic( H q ^ po g = _gQ( 5 ) x 1Q -11 (g) 
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and the light by light hadronic amplitudes @, |J 

a* adronic (light by light) = 8(9) x KT 11 



Altogether, one finds 



„ Hadronic 



6942(153) x 10 



-li 



(7) 



(8) 



Now we come to the electroweak contributions to a„, the primary focus of this 



work. At the one loop level, they are well known [10], [II], 12], [13], [L4 
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(9) 



where = 1.16639(1) x 10~ 5 GeV" 2 , M = M w or M ffigg s, and the weak mixing angle 
is defined by sin 2 6 W = = 1 - M^/Mf. We can safely neglect the O (m 2 /M 2 ) 
terms in (|9]). Also, throughout this paper we neglect terms suppressed by the factor 
1 — As^ whenever it simplifies the expressions without affecting accuracy. 

The one loop estimate of electroweak effects is about 5 times the anticipated 
experimental accuracy. Naively, one would expect higher order electroweak contribu- 
tions to be of relative 0(a/7r) and hence insignificant. However, an interesting study 
by Kukhto, Kuraev, Schiller, and Silagadze fl5| (KKSS) found that not to be the 
case. They showed that 2 loop electroweak contributions are quite large and must be 
included in any serious estimate of a^ w or confrontation with future experiments. 

Including 2 loops and making the approximations mentioned above, a^ w becomes 
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(10) 



It is natural to separate the subset of the two loop electroweak contributions which 
contain a fermion loop 



c = C lc " a + c 
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Parts of both (7 fcrm and C bos have been calculated by KKSS. Denoting the non 
calculated contributions by Rf (for fermionic loops) and by Rb (for the remaining 
diagrams) the KKSS results can be written as 



c 



fcrm 



c 



bos 





+ R b 



(12) 



The known parts reduce the electroweak contribution by about 24% (—46 x 
10~ n ), a significant decrease. A full calculation of Rb is quite a daunting task because 
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of the large number of diagrams. It has been estimated by KKSS to influence the 
correction factor C at the level of 10%. However, only a full 2 loop calculation will 
tell us if that is the case. 

In the present paper we reexamine the fermionic loops contributing to the two 
loop electroweak corrections and calculate Rf. We find that a significant subset of 
diagrams has been neglected in previous studies. In particular we find that the large 
logs of the ratios of Mz and lepton masses contributing to C ferm are cancelled by the 
corresponding quark diagrams. We also obtain new relatively large non-logarithmic 
corrections of 0(G^am^) and 0(1) terms. 

In our calculation we chose the ratio of muon and vector boson masses as an 
expansion parameter in the calculation of diagrams contributing to g — 2. Such 



asymptotic expansions have recently obtained firm theoretical foundation [16] . After 



the expansion we still have to perform two loop integrals, which however contain at 
most one mass scale. The calculation of such integrals is further facilitated by the 
integration by parts method [I7| and the symbolic manipulation programs written 
in FORM 0]. In some cases we used packages SHELL2 |L9| and MINCER |I| to 



check our results. 

We use dimensional regularization with the dimension of space-time equal D = 
4 — 2e and neglect terms containing 7^ and ln47r/i which accompany the poles 1/e 
and vanish in the sum; this explains the appearance of logs of dimensionful quantities 
in the intermediate (divergent) results. 

For the discussion of hadronic loops we use the following quark masses m u = = 
0.3 GeV, m s = 0.5 GeV, m c = 1.5 GeV, m b = 4.5 GeV, m t = 176 GeV. We perform 
the computations in the 't Hooft-Feynman gauge. The basic two loop diagrams with 
fermion loop contributing to muon g — 2 are shown in fig. 1. In addition we have to 
consider diagrams obtained by replacing vector boson propagators by corresponding 
Goldstone bosons as well as mirror counterparts of asymmetric diagrams. 

Contributions of diagrams with a fermion loop connected to the muon line via 
two charged bosons are shown in fig. l(a-c). Isospin +1/2 fermions are denoted by 
u and the isospin -1/2 fermions by d. We first consider the case when fermions in 
the loop belong to the first two generations. Here the masses of the fermions in the 
loop do not influence the result very much and we neglect them. The ratio of the 
neglected terms to the result is at most of the order of (ml/M^) ln(m^/ 'M^) < 0.3%, 
with m c denoting the mass of the charm quark w 1.5 GeV. In ref. [fJI| it has been 
argued that the diagrams 1(a) and 1(b) vanish by virtue of Furry's theorem. We find 
that this is not true even after adding contributions of all fermions in a generation. 
Furry's theorem consists in the observation that contributions of diagrams with two 
different orientations of the fermion loop mutually cancel. This is not the case for 
diagrams 1(a) and 1(b) because for every fermion flavor interacting with the external 
photon there is only one possible orientation of the fermion line. Only those parts of 
the expressions which contain a single 75 cancel out after adding contributions from 
the up-type quark, down-type quark and from the lepton. 

Adding contributions of quarks and leptons we obtain for a single light generation 
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(we put the relevant Kobayashi-Maskawa matrix elements equal 1) 



A /^rfcrm 
^ ° light 



(13) 



In the third generation we can neglect only masses of r and of the b quark. For 
the r lepton loop we obtain 
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(14) 



For the sum of all diagrams containing top and bottom quark loops we find 
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We computed 0{My V /m 2 ) and 0{M^/mf) corrections to this formula. Those terms 
turn out to have small coefficients which render them numerically insignificant. The 
singular terms m 2 /e will be canceled by renormalization of the W boson mass present 
in the one loop electroweak contributions to muon g — 2. 

We now consider the diagrams with a photon-photon- Z coupling induced by 
fermion loops shown in fig. 1(d). For electron or muon in the loop we obtain 



ACfr(e) 




(16) 



in agreement with ||15|| . In that reference one also finds a formula for the r lepton 
which can be generalized for all fermions sufficiently heavier than the muon and 
lighter than the Z boson 
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In practice this formula can be used for all quark loops except for the top quark, for 
which we find 



ACjT(f) 
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In ref. 15 the total fermionic two loop effect on muon g—2 was estimated by summing 



only electron, muon, and r contributions to diagram 1(d). It has been concluded that 
the source of large corrections are logarithms of ratios of these light fermion masses 
to the mass of the Z boson. Such treatment is incomplete and misleading. We can 
see from the formula (1171) that the sum over all fermions (see discussion in [HI) in 
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the first two generations leads to the cancellation of Mz dependent logarithms, due 
to the no-anomaly condition J2f^3fQ} = 0- This pattern no longer holds for the 



third generation; here, due to the large mass of the top quark, its contribution is 
suppressed by a factor M z jm\. This leads to the appearance of the logarithm of Z 
boson mass in the sum of all contributions to 1(d) 
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The first line of ( |19D gives the dominant contributions of the diagram 1(d) from all 
fermions. We note that the mass of the top quark is absent in this part reflecting 
the suppression of top loop discussed above. The second line summarizes the cor- 
rections to the dominant effect from the electron, muon, and from the top quark. 
Numerically we obtain Aai^ = —14.4 x 10~ n in contrast with the value given in 
21, p~5|| Aa lrf ( e MjT ) —25.6 x 1CT 11 . This reduction of the correction is caused by the 
cancellation of the Mz dependent logarithms. We stress that the Mz which is still 
present in the main part of ([19]) is caused by the large mass of the top quark and 



suppression of its contribution. 

To summarize this part we note that the large numerical value of the sum of 
diagrams 1(d) is generated by large mass splittings among the fermions. This is the 
main difference between our result and the result of [15|| : imagine a model in which all 
fermions had equal masses, then the sum of the three leptonic contributions discussed 
in |T5| would be equal 3 x Aai^), whereas we find that the total correction due to 
1(d) would vanish. We believe to have found a qualitatively new type of the top quark 
effect: it is namely the extremely large mass of the top quark which determines the 



shape of the main part of the formula fllTf), although the numerical value of m t is 
completely irrelevant. What is important is that M\jm? t <C 1. 

We now discuss the remaining diagrams, neglected in |I5| . In the diagram 1(e) 
we have to distinguish two cases, again treating top quark separately. For the light 
fermions we find, in contrast to the diagram 1(d), that the logarithmic factors are 
suppressed by extra powers of m*/M§; we retain good accuracy by taking massless 
fermions. We find 
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This contribution becomes sizable after adding the top quark effect. For all fermions 
together we find 
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where we neglected terms 0(M z /mf) which proved to be small. 
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In the remaining diagrams we have a scalar particle coupling to the fermion loop, 
and therefore we only consider top loops. If the Z boson in the diagram 1(d) is 
replaced by the neutral Goldstone boson we obtain 



^ A °ld(G) 



16 8 i m t 
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(22) 



Finally, there is the diagram 1(f) containing the Higgs boson. It is the only non- 
negligible contribution of the Higgs boson among two loop fermion diagrams. We 
consider three cases, depending on the hierarchy of masses of the top quark and the 
Higgs. For Mh -C m t we get 



AC{7 m 
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-2.1 • KT 11 (for M H = 60 GeV) (23) 
For the case of M H ^> m t we observe stronger suppression of this amplitude 
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(24) 



In order to estimate the size of the Higgs diagram in the case of similar top and Higgs 
masses we put m t = Mh- In this case we find 



ACf/ m (t) 
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(25) 



where CI2 is the Clausen function |22| . The contribution of the Higgs boson is small 
and we approximate it by — 1.5(±1.0) x 10~ n . 

A few words are in order to explain why several diagrams with neutral bosons 
have been omitted in fig. 1. Diagrams with two scalar bosons {HH, HG°, G°G°) 
coupling to the muon line are at most of the order m^/M|. The remaining diagrams 
(e.g. ZH) are either exactly zero, or are suppressed by the vector coupling of the Z 
boson to the muon (factor 1 — 4s^). 

Some large two loop corrections can be absorbed in the one loop result if it is 
parametrized in terms of determined from the muon's life time. This corresponds 
to the replacement of bare parameters 



8s 2 w M^ 



%1 
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(26) 



where A is determined by studying electroweak corrections to the muon decay width. 
Because we are interested in fermion loops, only quark and lepton corrections to 
the W propagator need be included. They induce the following counterterm which 
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cancels the divergences we encountered in the charged boson diagrams (eq. |Tj|) and 
the Z boson vacuum polarization (eq. pip 



ACc! m - + (1 2s^ + 

5 iy < j; 'wHy v ° 



3 m? / 1 1, 9 1, lf9 , 3 , , 79 \ 
— 5 1- h - In m, 2 - - In ME + In ML In m 2 

As^M^K 2e 2 * 5 z w 10 60 / 



The final result for the fermion loop effect on muon g — 2 is obtained by adding 
the contributions given by equations ( ]T3]JT^JT5|JT9||21| , p^ ,p^) and the contribution of 
the Higgs boson diagram ACfy m (t) given by approximating eq. p4|^5|) . Our final 
formula is 



5 (rndrn s rnb) 1 / 3 rn'f l rri T 

m? 3 , m? 

In ' 



16 S 2 W M^ 10S 2 W Myy 

A In U% - 11 _ _L_ + + AC{ c ; m (t). (28) 
5 Mf 5 lOs^ 15 x/ W V ; 

This is our main result which replaces the old estimate of C ' given in ([12]). We 
dropped all terms suppressed by negative powers of m t . We checked explicitly that 
their numerical impact is negligible. 

The 0(m 2 /M^) term in eq. (p8|) is related to the p parameter that appears 
in the ratio of weak neutral to charged current amplitudes and comparisons of the 
W ± and Z masses. It can be viewed as an induced correction brought about by 
our renormalization of the one loop Z contribution in terms of G^, a charged current 
parameter. We also note that except for their incomplete cancellation in the anomaly 
diagrams of fig. 1(d), no other effect of the 2 light fermion generations resides in our 
final result. 

For the numerical evaluation of the remaining terms which contain the weak 
mixing angle we take s 2 ^ = 0.223. We obtain 

which means that the correction to muon anomalous magnetic moment a M from the 
fermion loops is — (23±3) x 10~ n . The theoretical uncertainty has several sources: the 
unknown mass of the Higgs boson, uncertainty in the masses of the light quarks which 
parametrize the hadronic effects, and the large experimental error in the present value 
of mt- Finally, higher order three loop contributions remain unknown. Altogether we 
estimate these effects to yield an uncertainty at the level of 3 x 10~ n , more than an 
order of magnitude below the predicted experimental precision. 

Including the fermionic two loop corrections and partial two loop bosonic effects, 
we obtain the updated theoretical predictions 

aj w = (152(3) + 0A5R b ) x 10~ n 
af = (116591802(153) + 0.45/4) x 10~ n (30) 
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What remains is to compute Rb and lower the hadronic loop uncertainty. Work on 
both is in progress. 

After completing this calculation we learned about ref. [23| which contains an 
analysis of quark contributions to the diagram in fig. 1(d). For the light quarks their 
numerical result obtained using the chiral perturbation theory is the same as our 
evaluation using constituent quark masses. The large difference between their final 
evaluation of the fermionic loops and our result is due to diagrams of fig. l(a-c,e,f) 
which were not considered in |2l 
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